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ABSTRACT 

The saturation mechanism of Magneto- Rotational Instabihty (MRI) is exam- 
ined through analytical quasilinear theory and through nonlinear computation 
of a single mode in a rotating disk. We find that large-scale magnetic field is 
generated through the alpha effect (the correlated product of velocity and mag- 
netic field fluctuations) and causes the MRI mode to saturate. If the large-scale 
plasma flow is allowed to evolve, the mode can also saturate through its flow 
relaxation. In astrophysical plasmas, for which the flow cannot relax because of 
gravitational constraints, the mode saturates through field generation only. 

1. Introduction 

It is well known that coUisional hydrodynamic viscosity is too small to explain the 
inferred rate of angular momentum transport in accretion disks (Shakura & Sunyaev 
1973). One possible source of anomalous viscosity is turbulence resulting from the 
magneto- rotational instability, or MRI (Velikhov 1959; Balbus & Hawley 1991). In the 
linear phase, this instability derives its source of free energy from the radial gradient of the 
rotational flow in the presence of a weak magnetic fleld, i.e., the MRI requires the presence 
of a magnetic field, but is linearly stable if the field is too large. The properties of MRI 
turbulence, and its effect on angular momentum transport, depend in a fundamental way 
on the amplitude of the fluctuations in the non-linear state. We expect this amplitude to 
be related to the non-linear saturation mechanism of a single MRI mode. 

A common saturation mechanism for linear instabilities is the relaxation of the free 
energy source through quasi-linear processes. For the case of the MRI, this would imply 
modification of the radial profile of the rotational velocity, but this mechanism may not be 
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available in accretion disks and other astrophysical settings because of strong gravitational 
constraints that maintain a Keplerian profile (V^p oc r~^/^). Instead, the MRI may saturate 
(cease linear growth) if the non-linear evolution of the instability generates a mean 
component of the magnetic field that is sufficient to stabilize the mode. 

We remark that for most plasma instabilities, if the energy source is fixed and unvarying 
(and nonlinear coupling to other modes is ignored), the mode will grow without bound, 
until limited by dissipation or some physical dimension of the system. In this sense, the 
MRI might be atypical. 

In this paper we investigate the linear growth and non-linear saturation mechanism 
and amplitude of a single MRI mode in thick-disk geometry, i.e., ^ ~ 1. Our primary tool is 
linear and non- linear MHD computation in cylindrical (r, 0, z) geometry which is periodic in 
the azimuthal (0) and axial (z) directiona. We solve a model initial value problem in which 
the inner and outer radial boundaries are impermeable, perfectly conducting, concentric 
cylinders that can rotate independently at specified rates. They are coupled to the internal 
flow by hydrodynamic viscosity. The initial mean (azimuthally and axially averaged) profile 
V^{r) is Keplerian, and in most cases is assumed to be maintained for all times by external 
forces. Perturbations to the equilibrium that depend on (r, 0, z) are then introduced, and 
evolve dynamically according to the single fluid MHD equations. Linear growth rates 
and eigenfunctions are found by integrating the linearized, single fluid, visco-resistive 
MHD equations forward in time until an exponentially growing solution is obtained. The 
saturation mechanism and amplitude of the mode are determined by solving the non-linear 
MHD equations for a single mode, with azimuthal mode number m and axial mode number 
k, beginning from a small initial perturbation. At finite amplitude, the mode will interact 
with itself to modify the mean background state, which in turn alters both its temporal 
evolution and its radial dependence. In this sense it differs from a quasi-linear calculation. 
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in which only the effect of the linear eigenfunction on the background is accounted for. 

Our primary result is that, when the mean flow profile is maintained as Keplerian, the 
MRI can saturate due to the generation of a mean component of the magnetic field that is 
sufficiently large to linearly stabilize the instability. The mechanism for field generation is 
by means of a mean electric field <E> = — <VxB> =a < B > (where < ■ ■ ■ > denotes 
the azimuthal and axial average; and V and B are the velocity and magnetic fluctuations) 
that is produced by the nonlinear interaction of the finite amplitude mode with itself, and is 
parallel to the mean magnetic field. This is the well known a-effect, the correlated product 
of velocity and magnetic field fluctuations (Moffatt 1978). It is a fundamental dynamo 
mechanism. This suggests that the MRI may operate as a dynamo for the generation 
magnetic fields in accretion disks and other rotating systems. This is consistent with earlier 
results indicating that the magnetic energy increases during the onset of MRI turbulence 
(Hawley, Gammie & Balbus 1996). 

Finally, when the mean flow profile is allowed to evolve, a single MRI mode saturates by 
relaxation of the rotation profile, as expected. This is also accompanied by the generation 
of a mean magnetic field. Thus, the MRI dynamo mechanism (the alpha effect) may 
be robust. Our computational results are supported by an analytic theory in which the 
quasi-linear stresses and mean parallel electric field are calculated explicitly in terms of the 
linear eigenfunctions of the MRI. 

The paper is organized as follows. Previous work and the model are described in 
sections [2] and [3] respectively. Nonlinear evolution of a single axisymmetric MRI mode is 
presented in section HJ for fixed Keplerian flow (section 14.11) and for evolving flow (section 
14.21 ). Using nonlinear single mode computations, the alpha effect for non-axisymmetric 
MRI modes is examined in section [5l In section [6], the alpha effect for axisymmetric and 
non-axisymmetric MRI is obtained using quasilinear calculations. We summarize in section 
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2. Previous Work 

In early work, the linear theory of the MRI has been investigated carefully, and in the 
past decade nonlinear MHD simulations of the MRI have been performed (see, e. g., Balbus 
& Hawley 1998; Hawley & Balbus 1999 and the references therein), including evolution of 
an MRI-unstable system to a turbulent state. An interesting feature of the turbulent state 
is that small-scale magnetic fluctuations arise, representing an increase in magnetic energy 
(Hawley, Gammie & Balbus 1996). Turbulence generated by MRI in global disk simulations 
is one of the current MRI research topics (Armitage 1998; Machida, Hayashi & Matsumoto 
2000; Hawley 2000; Hawley 2001; Machida & Matsumoto 2003; Reynolds & Armitage 2001; 
Armitage & Reynolds 2003; McKinney & Gammie 2002; Nishikori, Machida & Matsumoto 
2006; Fromang & Nelson 2006; Lyra et al. 2008). 

Recently the issue of numerical convergence has been raised in the simulations of the 
MRI. It was found that the results of some simulations are not independent of the numerical 
dissipation scale, and physical dissipations should be used (Fromang & Papaloizou 2007). 
Numerical simulations of the MRI-induced turbulence in a shearing box with nonzero net 
vertical flux (Lesure & Longaretti 2007) and zero net flux (Fromang, Papaloizou, Lesur & 
Heinemann 2007) showed that the physical dissipation, the viscosity and the resistivity, 
affect the saturated level of MRI-driven MHD turbulence, and thus the amount of angular 
momentum transport. It was found that the turbulent activity is an increasing function of 
the magnetic Prandtl number. The dependence of the saturation level of the MRI on gas 
pressure has been investigated by Sano et al. 2004 using the local approximation. Based 
on series of local numerical simulations, a scaling law was also provided that describes the 
transport of angular momentum in shearing MHD boxes (Pessah, Chan & Psaltis 2007). 
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In addition to numerical simulations, analytical methods were used to study some 
nonlinear aspects of the MRI by Goodman & Xu 1994. Goodman & Xu introduced the 
notion of whether the MRI can amplify very weak initial fields until they reach equipartition 
with the gas pressure. They showed that MRI axisymmetric modes are exact nonlinear 
solutions of the local incompressible MHD equations, suggesting that the MRI perturbed 
magnetic field may grow to an amplitude many times larger than the initial unperturbed 
field. Their analytical calculations assumed a very special form for the axisymmetric MRI 
mode, with a purely vertical wavevector whose field perturbations are purely horizontal 
(perpendicular to the rotation axis). In this paper, we solve the full nonlinear MHD 
equations for a single MRI mode with perturbations in all three directions (radial, azimuthal 
and vertical). We find that large-scale mean magnetic field is amplified through the alpha 
effect that causes the mode to saturate. This effect was not studied by Goodman & Xu. 
In fact, the assumption made by Goodman & Xu may have limited the possibility of the 
MRI mode saturation through the alpha effect, which is the central result of our paper. 
In the limit of large amphfication of MRI perturbation, Goodman & Xu also showed that 
the MRI mode could themselves be subject to additional parasitic instabihties (such as the 
Kelvin- Helmholtz instability) which may destroy the primary MRI modes. The growth rate 
of parasitic instabilities is proportional to the MRI perturbation amplitude. Simulations of 
these instabilities would require multiple wavenumbers, and cannot be captured in single 
mode (single kz) computations. 

However, the saturation mechanism, a fundamental feature of MRI, is not fully 
understood. Most of the previous studies on the MRI saturation have dealt with what 
controls the saturated level of the MRI in a turbulent state. The most recent results show 
the effect of physical dissipation parameters on the saturation level using an unstratified 
shearing box model. However, these studies have not fundamentally investigated the 
cause of the mode saturation (Lesure & Longaretti 2007; Fromang, Papaloizou, Lesur & 
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Heinemann 2007; Pessah, Chan & Psaltis 2007; Sano et al. 2004). In this paper, we identify 
and explain the saturation mechanism of the MRI at the level of a single mode. Through 
the nonlinear evolution of a single mode, we simplify the problem to better identify the root 
of MRI saturation. Using this simple dynamical model, we find that the instability itself 
generates a large-scale magnetic field (mean magnetic field) through the alpha effect that 
causes the mode to saturate. 



3. The model 

We employ the MHD equations in doubly periodic (r, 4>, z) cylindrical geometry for 
both the analytic (quasi-linear) and computational studies. All variables are decomposed 
as fir,cl),z,t) = Eim,k)fmAr,t)e'^""^+'''^ =< fir,t) > +fir,(f) z,t), where < / > is the 
mean {m = k = 0) component, and / is the fluctuating component (i.e., all other Fourier 
components with m ^ and k ^ 0). Note that the mean component is a function of 
radius. We consider an azimuthal equilibrium flow Vq = V^(r)0 and vertical magnetic field 
B = Bqz. (Here, V^(r) is the equilibrium value of the mean azimuthal flow < >, and Bq 
is a constant.) The single fluid MHD equations are, 
8A 

^ -E^SVxB-rjJ (1) 

dt ' ^ ^ 

= -SpV.VV + S3 xB + P^V^V-S^VP-SpV^ + Fg (2) 

^ = -5v • (pv) - 5(r - i)pv ■ V (3) 

1^ = -SV-{pN) (4) 
B = VxA (5) 
J = VxB (6) 

where the variables, p, P,V,B, J, F, and <I> are the density, pressure, velocity, magnetic 
field, current, ratio of the specific heats, and gravitational potential, respectively. Time 
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and radius are normalized to the resistive diffusion time tr — 47rr2^/ 0^770 and the outer 
radius r2, making the normahzed outer radius unity. Velocity V and magnetic field B are 
normalized to the Alfven velocity V^, and the magnetic field on axis Bq, respectively The 
parameter 5" = — is the Lundquist number (where ta — ^2/!^), and = — measures 
the ratio of characteristic viscosity vq to resistivity r)o (the magnetic Prandtl number), 
where t^s = 47rr2^/c^t'o is the viscous difi'usion time. The factor (3q = 8itPq/Bq is the beta 
normalized to the axis value. The resistivity profile rj is uniform. The boundary conditions 
in the radial direction are as is appropriate to dissipative MHD with a perfectly conducting 
boundary: the tangential electric field, the normal component of the magnetic field, and the 
normal component of the velocity vanish, and the tangential component of the velocity is 
the rotational velocity of the wall. The azimuthal (0) and axial (z) directions are periodic. 

We pose an initial value problem that consists of the equilibrium plus a perturbation 
of the form /(r, 0, z,^) = /^^^(r, i)exp(im0 + ikz). Equations (1-6) are then integrated 
forward in time using the DEBS code. The DEBS code uses a finite difference method 
with a staggered grid for radial discretization and pseudospectral method for azimuthal 
and vertical coordinates. In this decomposition, each mode satisfies a separate equation of 
the form dfm,k/dt = Lm,kfm,k + Z)(m',fe') ^m,k,m',k'^ where Lm,k is a linear operator that 
depends on fofi{r,t), and A^^^^/ ^/ is a nonlinear term that represents the coupling of the 
mode {m,k) to all other modes {m\k'). (This latter term is evaluated pseudospectrally.) 
The time advance is a combination of the leapfrog and semi-implicit methods (Schnack et 
al. 1987). 

We distinguish between three types of initial value computations. In linear 
computations, the initial conditions consist of an equilibrium < /(r) > plus a single 



mode /m.fc(r, O)ea;p(im0 + ikz), as described above, with fm,k/ < f > 



<< 1. The initial 



amplitude fm,k{'>^, 0) is a polynomial in r that satisfies the boundary conditions at r = ri 
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and r — r2- The initial amplitudes of all other modes are set to zero, and are constrained 
to remain zero for all time. Only the mode (m, k) is then evolved; in particular, the 
equilibrium (the m = 0, k = mode, or < / >) is not evolved, and remains fixed in time. 
The mode (m. A;) then satisfies the linear equation dfm,k/dt — Lm,kfm,k, where Ljn,k{f) is 
independent of time; an exponentially growing solution indicates linear instability, and the 
solution fra,k{'^'-i't) IS the radial profile of the unstable eigenfunction. 

In nonlinear single mode computations, the problem is initialized in the same 
manner. However, the m = 0, /c = component (the background) is allowed to evolve 
self-consistently. The term A^o,o,m,ifc contains the nonlinear contribution from the mode 
(m, k) interacting with itself, and the operator which depends on /o,o, can evolve with 
time. All other modes (m', k') have and retain zero amplitude. When the mode amplitude is 



small 



fm,k/ < f > << 1)) the nonlinear term is negligible and the linear solution is again 



obtained. As the unstable mode grows and 



fm,k/ < f > 



becomes finite, the evolution of 



the background profile /o,o can modify the term Lm,k{f-it) and this can affect the evolution 
the mode (m, k). If, as a result, the mode (m, k) ceases to grow exponentially it is said to 
saturate. (In most cases, this means that the mode ceases continual growth altogether.) We 
remark that this differs from quasi-linear saturation, because in the saturated state the 
radial profile of the finite amplitude solution fm,k{''^,t) can differ significantly from that of 
the linearly unstable eigenfunction. 

In a fully nonlinear computation, all modes are intitialized with small random 
amplitude and are evolved in time, including the full nonlinear term km! k' ■ this case, 
saturation can occur as a result of both background modification and coupling to other 
modes. If there is sufficient numerical resolution, this can lead to a turbulent state. 

In this work we consider both linear and nonlinear single mode computations of the 
MRI. The latter begin from a state in which the plasma rotates azimuthally with a mean 
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Keplerian flow < V^(r) >oc r~^/^. The initial (i.e., at t = 0) radial equilibrium force balance 
(Eq. 2) is satisfied by ^Vp + pV$ = pV^/r, where V$ = GM/r'^, and the initial pressure 
and density profiles are assumed to be radially uniform. 

The azimuthal fiow profile is the source of free energy for the MRI. We might expect 
that the MRI would saturate by altering this profile. However, in astrophysical settings the 
Keplerian profile is fixed by gravitational constraints, so this saturation mechanism may not 
be available. We investigate the saturation process in this case by performing computations 
in which the mean fiow profile is fixed in time (see section H7T1) . In these computations, the 
mean Keplerian profile is maintained in time by an external force; an ad-hoc force (F^) is 
added in the momentum equation to fix the mean fiow profile. In other computations we 
allow the mean fiow profile to evolve in time (see section l42ll : the initial Keplerian fiow 
consistent with the angular velocity at the inner and outer boundaries is allowed to relax 
self-consistently as the mode reaches finite amplitude. This is the usual case of relaxation 
of the source of free energy for the mode. In both cases, all other mean variables (e.g., 
the mean magnetic field profile) are allowed to evolve. Pressure and density are evolved; 
however, they remain fairly uniform during the computations. (Note that perturbations to 
the azimuthal velocity are allowed to evolve.) 

In all cases the MRI is found to cease exponential growth, and the background profile 

< /(r) > is modified. The saturation mechanism for a single MRI mode is then investigated 
by performing linear computations in which the initial background state /o,o('";0) is 

the saturated background profile obtained from the nonlinear single mode calculations, 

< /(r) >. We find that this new linear solution does not grow exponentially, indicating 
that, as a result of the nonlinear single mode evolution, the background state has become 
linearly stable. In particular, this state has a modified mean magnetic field profile that has 
sufficient magnitude to stabilize the mode. 
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Fully nonlinear calculations of MRI saturation will be reported in a future publication. 

We consider a cylindrical disk-shaped plasma with aspect ratio L/(r2 — ri) (where L 
is the vertical height, and ri and r2 are the inner and outer radii, respectively). The inner 
and outer radial boundaries are perfectly conducting, concentric cylinders that can rotate 
independently at specified rates. Periodic boundary conditions are used in the vertical and 
azimuthal directions. 

The aspect ratio (L/ (r2 — Ti)) used in the nonlinear computations is 1.3. The nonlinear 
computations are performed in a thick-disk approximation where vertical and radial 
distances are of the same order. The range of parameters used in the computations is 

/5 = 10^ - lo^ s = 10"^- lo^ Pm = 0.1 - 1. 

4. Axisymmetric single mode nonlinear evolution 

We present the results of single mode nonlinear computations with both fixed and 
evolving mean azimuthal velocity profiles for an axisymmetric MRI mode (i.e., m — 0), 
as these modes have the largest linear growth rate (see, e.g., Balbus & Hawley 1998). 
The fixed mean flow profile is relevant to the astrophysical setting; evolving mean fiow is 
the usual case of relaxation of the free energy source. In both cases the initial condition 
consists of the equilibrium plus a perturbation with single vales of m (in this case m — 0) 
and k. In the first, we maintain a fixed mean profile < >, so as to isolate the effect of 
magnetic field growth on mode saturation. In the second, the mean azimuthal fiow profile 
is permitted to evolve in time. (In both cases, perturbations to the azimuthal fiow are 
permitted and evolve.) We find that mean magnetic field growth plays an important role 
in the saturation of the mode, although the effect is somewhat stronger when the mean 
azimuthal flow is flxed in time. 



4.1. Evolution with fixed mean flow 



We computationally solve the MHD equations (section [3]) with initial conditions of 
Keplerian flow (V^ = Vor~^/^) and a weak vertical magnetic fleld B = Bqz. We flrst, 
consider only axisymmetric modes (m = 0), which are the fastest growing modes and will 
be treated analytically in section [61 Since the computations are single mode (single m and 
k), the spatial resolution is only in the radial direction. The number of radial grid points 
used is 250. The solutions are converged in timestep and radial spatial resolution. We have 
performed nonlinear computations with different axial wave number and different Lundquist 
numbers, 5* = 10^, 10'^, 10^, and with different equilibrium azimuthal flow magnitudes 
Vo = 10 — 100 Va- We flnd that the dynamics of mode saturation (and fleld ampliflcation) 
is independent of axial wave number k and flow amplitudes, and very similar for the three 
different Lundquist numbers. Therefore, we only present nonlinear computation for the 
axisymmetric k = 6 mode with Vq = 16Va, at S = 10^ and S = 10^. 

Figure [1] shows the components of the magnetic modal energy during the nonlinear 
evolution when flow is flxed for S* = 10^ and S = 10^. The axisymmetric mode grows 
with the growth rate jta = 9.5 and starts to saturate around t=0.016 (Fig. [11(a)). Kinetic 
energy at t=0 is about two order of magnitude (50 times) larger than magnetic energy; 
however after saturation of the mode, total magnetic energy increases, and it becomes of 
the same order as the kinetic energy. Mean vertical magnetic fleld at t=0 is = 1, and is 
shown with the dashed-dotted line in Fig. [21(a). As the mode starts to saturate, the vertical 
magnetic fleld starts to increase as the result of the alpha effect < V x B >^. Figure [21(a) 
shows the mean vertical magnetic fleld at time t=0.016 for S = 10^ (around the time of 
mode saturation). At this time the vertical mean magnetic fleld in the core (r=0.2-0.4) is 
increased by a factor of nine. The azimuthal component of the alpha effect which causes the 
generation of the vertical magnetic field is large in the core (Fig. [21(b)). The mean vertical 
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magnetic field and the azimuthal component of the alpha effect are also shown in Fig. [2] for 
S = 10^ (dashed line) around the time of mode saturation. Similar radial structures and 
magnetic field amplification are observed for S = 10^. The radial structure of the nonlinear 
alpha effect (Fig. [2](b)) is consistent with the quasilinear calculations (section 6, Fig. [1]) and 
results in the amplification of magnetic field. Therefore, the nonlinear computation shows 
that the fiuct nation- induced alpha effect amplifies the mean magnetic field in the core and 
causes the saturation of the mode. 

We verify the stabilizing effect of the modified mean magnetic field profiles by 
performing linear computations of the saturated state, as described in section 3. These 
computations use the background profiles near the time of saturation. In this way we 
are able to calculate the linear growth rates of any resulting instabilities. Each linear 
computation is initialized with perturbation of a single Fourier mode (m=0 and k=6) 
and the equilibrium profiles. For each linear computation, the equilibrium profiles are the 
mean (m=0 and k=0 component) magnetic and velocity fields taken from the single mode 
nonlinear computation at a specific time near the saturated state (e. g. ti marked in 
Fig.EKa)). 

We perform four separate linear computations for m=0, k=6 mode which are initialized 
with four different equilibrium profiles. Four set of equilibrium profiles are the mean fields 
taken from the single mode nonlinear computation at four times during the saturation. 
Figure [3](a) shows the radial magnetic energy vs. time for the nonlinear computation with 
fixed mean fiow; the four times at which the linear computations are performed are also 
marked. The linear growth rates calculated from the four linear computations decrease as 
the mode saturates from ti to t4. The linear growth rate decreases from 0.58 (the value 
before the saturation at t=to) to 0.048 at t4. The mean vertical magnetic field at these 
times is also shown in Fig. [3]^b). As the mean vertical magnetic field amplifies from ti to t4 
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by the alpha effect in the plasma core, the resulting mode growth rate obtained from the 
linear computation substantially decreases (Fig. 3 (c)). 

We find that the dynamics of the saturation and the resulting alpha effect is 
independent of axial wave number k. To illustrate this effect, we have performed nonlinear 
computation of a single mode with higher k (k=30) shown in Fig. 3 (d). The mode with 
k=30 grows faster than k=6 with a growth rate '^Torut = 1-17 and starts to saturate 
around t=0.007. The saturation level for this mode is lower than k=6 mode. Similarly, 
we initialize and perform linear computations with the saturated profiles from nonlinear 
computations at times ti - t^. The saturated mean vertical magnetic fields which are used 
as the equilibrium profiles for the linear computations are also shown in Fig. 3(e). Three 
separate linear computations are initialized at times ti-ts shown in Fig. 3(d). Because the 
linear computations have been started with the saturated mean fields, the linear growth 
rates obtained from these three computations are zero or negative, as seen in Fig. 3(d). 
The linear stability computations thus show that the amplification of the mean vertical 
magnetic field by the azimuthal component of the alpha effect < V x B >^ suppresses the 
mode growth and thus causes the mode saturation. The a-effect of the saturated MRI 
therefore appears to be independent of the axial mode number k. 



4.2. Evolution with evolving mean flow 

In the computation with a fixed mean azimuthal flow profile (Section 4.1), the mean 
magnetic field profile could be modified by the fluctuation-induced alpha effect. In the 
case when mean flow is allowed to evolve, the mean azimuthal flow profile will also be 
modified by the MHD stresses. Figure H] shows the nonlinear evolution of m=0, k=6 mode 
for evolving flow. Similar to the computation with fixed flow (Fig. [1]), axisymmetric mode 
grows and starts to saturate around t=0.016. Figure [5] shows the evolution of magnetic field 
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and flow during the nonlinear computations at times t = to (initial), ^1,^2 ( near saturation) 
and ta (final saturated). As is seen, the mean flow gradient is reduced in the core during 
the saturated state. The initial Keplerian mean flow profile flattens in the core (Fig. El (b)). 
Thus, with smaller mean flow gradient (less free energy for the mode growth), the mode 
starts to saturate. In addition to the mean flow modification, magnetic field generation still 
occurs. The magnetic field amplifies by a factor five in the core (FiglSJ^c)) as the azimuthal 
component of the alpha effect (Figj5](d)) increases. Since the mean flow is allowed to evolve, 
the amplification of the magnetic field is smaller than the computation with fixed mean 
flow. 

The alteration in mean flow and mean magnetic field profiles both contribute to mode 
saturation. The Reynolds and Maxwell stresses responsible for the momentum transport 
that flattens the flow profile are shown in Fig. [21(a). Both stress terms contribute to the flow 
relaxation and outward momentum transport (consistent with Eqs. [18], [22] in section 6). The 
ratio of the Maxwell and Reynolds stresses depends upon the magnetic Prandtl number. 
Fig. Mjo) shows the profiles of Maxwell and Reynolds stresses for the same computations 
but with Pm=0.1 (final saturated state). It is seen that the Reynolds stress is larger than 
the Maxwell stress in the core for computation with lower Prandtl number. Final relaxed 
flow profiles for these computation with Pm=l and Pm=0.1 is slightly different in the core 

(Fig. [ZD- 

5. Non-axisymmetric single mode nonlinear evolution 

Finally, we investigate the generation of magnetic field by the alpha effect from 
non-axisymmetric MRI modes. Nonlinear computation for an MRI mode with m = 1 is 
performed with fixed Keplerian mean flow profile with the same equilibrium and parameters 
as the m = mode. The non-axisymmetric m=l mode grows more slowly than the m = 
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mode, but saturates at the same level (Fig. [8]). We find that the saturation mechanism 
for the m=l mode is very similar to that of the m=0 mode. However for the m = 1 
mode, the vertical component of the alpha effect, < V x B >z, is also nonzero and 
becomes important in the process of mode saturation. The vertical component of the alpha 
effect will be calculated in the next section (Eq. [12], [T3l) and shown to be nonzero for the 
non-axisymmetric MRI modes. Figure MJo) shows the radial structure of both vertical and 
azimuthal components of the alpha effect for the m=l mode during the linear growth. 
Because of the nonzero vertical alpha effect, an azimuthal magnetic field i?,^ (which is zero 
at t=0) is generated and is shown in Fig. [9l Thus, the saturation of m=l mode occurs 
through amplification of both vertical and azimuthal magnetic field through the alpha 
effect. 



To evaluate the infiuence of the instability (fiuctuations) on the mean magnetic field, we 
calculate the alpha effect — the mean fluctuation-induced electromotive force < V x B >, 
where V and B are the velocity and magnetic fluctuations, and <> denotes an average 
over the azimuthal and vertical directions. Perturbed quantities are assumed in the form of 
exp{'-ft + im(f) + ikz). The alpha effect is obtained for a single MRI Fourier mode (single m 
and k). The linearized momentum equations in terms of the Lagrangian displacement ^ are 
(Chandrasekhar 1961), 



6. Quasilinear calculations 



_dX_ 

dr 
imX 



(7) 



8 



r 



(9) 
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where 7 = 7 + imO,, X = P + Bq ■ B, lua = k B^/p and Q, = V^/r is angular velocity. 
The velocity fluctuations are related to the Lagrangian displacement ^, K = 7^r5 

= 7^(^ — — -y-J^r and Vz = 7^^. In an ideal plasma (zero resistivity) magnetic 
field perturbations are also directly related to the displacement, B = ikB^^. Using 
incompressibility (^^ = —l/ikr{r-^r)' — m/rk^^) and Eqs. 019], azimuthal and vertical 
displacements are obtained in terms of C,r, 

1 ' -'"^ & + (10) 



'1 + m2/(r2/c2) (^2 _|_ 1^^-) ^ ^2]^2 

^ _ —mkr . —2^7 ^ _|_ ^ ^^^^ -^^^ ^ yi ^^^j^-j 

^ (r^/c^ + m?) (7^ + o;^) *" rk [r'^k'^ + m^) ^' 

We then construct the quasilinear alpha effect , <V x B >z= l/2Re(y*B^ — V^Br) and 
<V X B >^= l/2Re{V*Br — V*Bz), in terms of the radial displacement ^r, 

< VxB >^= ^Ur^r)'- , ,M rir)'+ ^(^i-f)C' (12) 



~ ~ m'jkB ^ , ^ 2m'ykBVt^ . ^ 2rrr''ykBVL'^ 
< ^ X i? >.= -^-^^^e.(re)' + (^i - l')e jf, C (13) 

where, 

H' = {1 + |^)(4mV^^' + (7' - m'^' + ^1)') (14) 

For axisymmetric m=0 modes, the alpha effect reduces to a simple form. The vertical 
component of the alpha effect (Eq. [T^ vanishes. The azimuthal component (Eq. [T^ is 
nonzero and becomes, 



<VxB>^= ^Vr-^{rVr) (15) 
7r or 

Thus, the instability is accompanied by an alpha effect that generates magnetic field. 

To obtain the radial structure of the alpha effect, linear computation (Eqs. 1-6) is 
performed for axisymmetric modes. We construct the alpha effect using linear eigenfunctions 
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obtained from single mode computation. Figure [TOT a) (solid line) shows the azimuthal 
component of the alpha effect for the m=0 mode with k = 6. The azimuthal component 
generates axial magnetic field in the inner region of the disk, a stabilizing influence. At 
higher k values (Fig. [TO](a). dashed line) the mode structure and alpha effect becomes 
more localized radially. Linear computations are also performed for different Lundquist 
numbers, S = 10^, 10^, 10^ (Fig- fTUT b)) . The growth rate of the MRI mode does not 
change significantly with Lundquist number (since MRI is an ideal mode). It is seen that 
the structure is slightly broader for lower S, which is expected because of higher dissipation 
at lower S. Similar radial mode structure (Vr) with peaked structure near the inner radius 
has also been obtained in earlier work by Curry, Pudritz & Sutherland 1994 and Blokland 
et al. 2005. 

We can similarly evaluate the fluid stresses that alter the flow. The MHD stresses for 
axisymmetric MRI mode are given by Balbus and Hawley 2002. Here we also calculate the 
quasilinear stresses for a non-axisymmetric MRI mode. The nonlinear evolution of this 
mode was examined in section [5l Using eqsOIHl the Maxwell stresses (for both axisymmetric 
and non-axisymmetric modes) are calculated in terms of radial displacement 

< r BrB^ >= — — (Ibj 



< rBrB^ >= — H — (17) 



It is seen that the vertical Maxwell stress is zero for the axisymmetric modes and the 
azimuthal Lorentz force < J x B >^= < f^B^B^ > is obtained 

<JxB>^=^^{DQ) (18) 
ar 

7(72 + a;i) 



Similarly, the Reynolds stresses are 



< >= -^2!5(^±^g _ ,rHiil - - '^l^C (20) 
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Again, the only nonzero component for m=0 modes is azimuthal, and < V ■ W >^= 
2^ J: < r^VrV^ > is calculated to be, 



< ^ ■ >*= ^^^m) + J^^ym-') (22) 

where the coefficient D is defined in Eq. [191 Adding eqns [18] and [22l the total azimuthal 
fiuid force for m=0 modes becomes, 

<JxB>,+p<V.VV>,= {e + ^,)^^{Dn) + J^^^y^'\Vr\') (23) 

The fiuid force consists of a simple diffusion-like part (the first term) and a more complex 
part (the second term) that also depends on fiow nonuniformity (fi'). 



7. Summary 

In summary, we have shown that the MRI will generate large-scale magnetic field as 
a mechanism for saturation. Computation of the fully evolved, nonlinear steady state of a 
single mode establishes that the mechanism that determines mode saturation is as described 
by quasilinear theory, and that magnetic field generation persists. This is an alpha effect, 
and can be understood through the stabilizing effect of the magnetic field. The analytic 
and computational work presented here treats one mode, as a means to clarify the effect. 
The influence of multiple modes and mode coupling will be of considerable importance. 
In a future work we will investigated the MRI saturation via the alpha effect using the 
full nonlinear computations with multiple modes. The fleld generated is large-scale in that 
it is symmetric in toroidal and axial directions. However, it varies radially such that the 
total toroidal and axial magnetic flux within the plasma does not change signiflcantly in 
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time. Thus, the examples studied here would not be defined as a "large-scale dynamo." 
Future work also will investigate whether there are geometric variations for which the 
generation of large-scale field would be accompanied by the generation of total magnetic fiux. 

We acknowledge useful discussions with Ellen Zweibel. We also thank the referee for 
useful comments. This work is supported by National Science Foundation. 
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Fig. 1. — Nonlinear evolution of axisymmetric MRI mode. Magnetic energies, 
= 1/2 / B^(m,nydr^ (solid line), Wr = 1/2 / Br.(m,n?dr^ (dashed line), and W, = 
1/2 J Bz(m,nydr'^ (dashcd-dotted line) vs. time for computation with mean flow profile fixed 
(/? = 10^ Pm^l)- (a) S = 10^ (b) S = 10^ 
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Fig. 2. — Radial profiles when mode starts to saturate in computation with fixed mean fiow 
of (a) mean vertical magnetic field (b) azimuthal component of the alpha effect, for S — 10^ 
(solid fine) at t=0.016 and S = 10^ (dashed fine) at t=0.0015 . 
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Fig. 3. — (a)-(c) Computations for m=0, k=6 mode, (d)-(e) Computations for m=0, k=30. 
(a) Radial magnetic energy during the nonlinear evolution of MRI mode (m=0, k=6) with 
fixed mean flow. Four times during the saturation are marked, (b) Mean vertical magnetic 
field profiles for these four times during the nonlinear saturation. Four linear computations 
are initialized with these mean magnetic field profiles, (c) Growth rates ( jTorbu) calculated 
from the linear m=0, k=6 computations, (d) Radial magnetic energy during the nonlinear 
evolution of MRI mode (m=0, k=30) with fixed mean flow. Three linear computations (at 
^i^^s) show no growth for k=30 mode, (e) Mean vertical magnetic field profiles for ti-t^ 
times during the nonlinear saturation used in the linear computations (m=0, k=30). 
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Fig. 4. — Nonlinear evolution of axisymmetric MRI mode. Total magnetic energy Wm,n — 
1/2 J B?^dr^ vs time for computation when flow profile evolves (/3 = 10^, Pm — 1, S — 10^). 
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Fig. 5. — Nonlinear m=0 mode computation with evolving flow, (a) averaged radial magnetic 
field fluctuations vs time, (b)-(d) radial profiles of azimuthal flow, vertical magnetic field 
and the alpha effect {S <V x B >^) at 4 times. 
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Fig. 7. — Final relaxed flow proflles from nonlinear computation for two Prandtl numbers, 
Pm=l and Pm=0.1. 
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Fig. 8. — Nonlinear evolution of non-axisymmetric m=l MRI mode, (a) Magnetic energy 
Wjn,n = 1/2/ -Sr(mn^'"^ time for Computation when mean flow profile is fixed for m=l 
(solid line) and m=0 (dashed hne) (/3 — 10^, = 1). (b) Azimuthal and vertical components 
of the alpha effect during the linear growth of m=l mode. 
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Fig. 9. — Radial profiles of vertical and azimuthal magnetic field for m=l mode at the time 
of saturation. Dashed line denotes mean magnetic profile (5^) at t=0. 
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Fig. 10. — Radial structure of the alpha effect < V x B >^ ior the m=0 MR mode (a) for 
two wave numbers k = 6, 60 (5 = 10^), (b) for S = 10^, 10^ 10^ with k = 6. V^/Va = 51, 
/? = 10^, Pm = 1. The growth rates are ■jTorbit = 0.23 and ■jTorbu = 1-7 for A; = 6 and A; = 60 
respectively {Torbu — Stt/JI). 



